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Abstract 

We describe the decomposition of one-loop QCD amplitudes in terms of colour-ordered building 
blocks. We give new expressions for the coefficients of QCD colour structures in terms of ordered 
objects called primitive amplitudes, for processes with up to seven partons. These results are needed 
in computations of high-multiplicity scattering cross sections in next-to-leading-order (NLO) QCD. 
We explain the origin of new relations between multi-quark primitive amplitudes which can be 
used to optimise efficiency of NLO computations. As a first application we compute the full-colour 
virtual contribution to the cross section for W -\- 4-jet production at the Large Hadron Collider, 
and verify that it is very well approximated by keeping only the leading terms in an expansion 
around the formal limit of a large number of colours. 
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I. INTRODUCTION 



Scattering processes at the Large Hadron Collider (LHC) are dominated by the strong 
interactions. Of particular interest in the search for new physics are those which occur at 
high momentum transfer. These processes frequently mimic new physics signatures, and so 
a solid theoretical understanding is important to maximize the sensitivity of experimental 
searches. Perturbative calculations in QCD, including quantum corrections, provide first- 
principle predictions in this direction, strongly linking fundamental theory to experiment. 

Recent years have seen remarkable progress in perturbative QCD, with new predictions 
at next-to-leading order (NLO) in the QCD coupling for many processes of interest at the 
LHC The present state-of-the-art studies of parton-level NLO QCD processes include 
three 2|-[6| or four 7H13| coloured final states. Such computations are challenging and 
rely on efficient algorithms for organising the various dynamic degrees of freedom, such as 
kinematics, spin and colour quantum numbers. In this article we focus on one particular 
aspect of these algorithms, namely the handling of the colour degrees of freedom for loop 
computations. Our results are of immediate use in modern approaches to NLO computations 
in QCD. 

At tree level various approaches for dealing with colour have been developed with differing 
advantages. The colour-ordered approach groups together terms with identical SU{Nc) 



or tree-level 
. This 



16 



colour factors (where is the number of colours). It was first developed 
amplitudes in IJ, |l5| and later extended for use in loop computations in refs. 
approach allows the symmetry properties and simplified kinematic dependence of ordered 
amplitudes to be exploited. Alternatively, one can treat colour simultaneously with the 



kinematic variables 



20|. 



Berends-Giele recursion 



'his colour-dressed method, most effectively implemented as a 



2l|, has been shown to be very efficient 



and is naturally applied to generic Standard Model processes 



br tree-level computation 



22|-|25|. 



At loop level the interaction of colour degrees of freedom is more complex, and requires 



refined techniques. For a large number o: 



coloured final states the ordered approach has 



proved particularly useful (see e.g. |3l-l5l.llll.ll2|). For this colour organization an expansion in 
powers oil/Nc is naturally set up which allows for significant efficiency gains in the numerical 
phase space integration [5|]; for a fixed integration error, fewer evaluations are needed for 
the smaller, but typically more time consuming, subleading-colour contributions 26 1. 
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The colour-ordered approach is also naturally suited to modern methods for loop compu 



tations using 



unitarity 
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unitarity 
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3, Si and on-shell recursio n | 29l-l3ll| . In particular, generalised 
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37| was first developed in the 



33| and its numerical extension 
context of the colour-ordered approach. Ordered loop amplitudes have simplified analytic 
properties and can be built through unitarity cuts from ordered tree amplitudes. The di- 
vision of the full amplitude into ordered pieces also leads to improved numerical control. 
By now a number of different methods for dealing with colour are in use for computing 



hard virtual scattering matrix elements 
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approaches 



IG 



43| . This includes colour sampling recursive 



391], in analogy to those developed at tree level. More traditional methods to 



organise loop computations usin g co lour information in the context of a Feynman- diagram 



calculation have been applied in 



44|. 



Here we will follow the method based on colour-ordered amplitudes. In this approach 
the colour information is factorised from g aug e- invariant purely kinematic pieces, the latter 
being embodied in primitive amplitudes [18!]. Intermediate steps of the computation are 
then free of any colour information. This information is processed separately and reincor- 
porated at the end of the calculation. As yet this approach is the only one which has led to 



phenomenological applications with five final state objects [11|, Il2|. However, one technical 
issue with the ordered approach is that the decomposition of matrix elements, and in partic- 
ular the virtual piece, in terms of primitive amplitudes, can be non-trivial, especially as the 
number of fermion lines increases. The central difficulty is that the fiow of quarks through 
a colour-ordered diagram is not uniquely determined by the external colour charges; each 
quark may navigate the loop in two distinct ways. In addition, the 1/Nc terms in the gluon 
propagator do not cancel when exchanged between quark lines, as they do in the purely 
gluonic case. In this paper we show how the derivation of the colour decomposition can be 
automated using Feynman diagrams. This only needs to be done once for a given process, 
and is therefore not an issue for the numerical efficiency of an NLO program. 

At leading order in the number of colours, the decomposition in terms of primitives is in 
general straightforward to write down. The corresponding results at subleading colour are 
non-trivial to derive, in particular, when multip le quar k lines are present; up to now they 
have been dealt with on a case-by-case basis 
will derive explicit expressions at fixed multiplicity for QCD amplitudes with many quark 
lines in terms of ordered objects. In particular, we present new results for the decomposition 
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4 71 ]. Using our automated approach we 
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of four- and six-quark amplitudes with up to seven coloured states. (See ancillary text files 
for our explicit expressions jsGj.) Thus we provide the missing pieces for the complete colour 
decomposition of QCD amplitudes with six and seven coloured states. These new results 
are of much wider relevance, however. With the colour decomposition of the pure QCD 
amplitudes at hand one can obtain rather simply their generalisations to include additional 
colourless objects such as leptons, vector bosons and Higgs bosons. 

As a state-of-the-art application of our results we present distributions of the virtual 



contribution to W + 4-jet production including su 



Dleading-colour terms. A study at the 



LHC of this process has been published in ref. [llj. In that study the finite part of the 
virtual piece was given at leading order in a large- A^c colour expansion. Here we confirm 
that the contribution of subleading-colour terms is small (a few percent), and shifts the 
virtual contribution uniformly over phase space. 

Finally, we discuss an additional result, that is certain relations among the set of multi- 
quark primitive amplitudes. These relations appear as a by-product of our method and 
originate in the antisymmetry of the colour stripped fermion-fermion-gluon vertex. While 
manifest at low multiplicity, these "fermion-flip" relations become more and more intricate 
as coloured states are added. From a practical point of view, such relations enable us to 
express one-loop amplitudes in different ways. This can be exploited to obtain a minimal set 
of contributing ordered amplitudes, allowing a reduction of the computing resources needed 
for NLO calculations. 

The organisation of this article is as follows. We begin in section [III by describing in detail 
our conventions for dealing with the colour degrees of freedom. In section IIIII we describe 
our algorithm and some details of our calculation. We also discuss how relations between 
primitive amplitudes are obtained as a by-product. In section |IV] we present selected results 
and a numerical study of the virtual contribution to 14^ + 4-jet production. Finally, we draw 
our conclusions and present an outlook for future research. 

II. COLOUR DECOMPOSITION 

In this section we review how full QCD one-loop amplitudes can be decomposed into 
contributions associated with SU(Nc) colour structures. While some of the notation extends 

nn 

the current literature, most notation and concepts follow the detailed discussions in p^, 118|, 
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19| and the reviews 48|, |49|. (See also ref. 50| for a recent review with topics related to this 
subject.) 



A. Basic conventions 

We will be concerned with QCD scattering amplitudes of multiple quarks and gluons. 
Quarks are denoted by q and q and are assumed to transform in the fundamental N^. and 
anti-fundamental Nc representations of the gauge group SU{Nc). When multiple quark 
flavours appear^, they will be distinguished by a flavour index qi and with I = 1,2, . . . Uf. 
In the absence of this index, the quarks q and q are understood to have a flavour number 
of one. We will keep the rank A^^ as a free parameter and specialize to A^c = 3 only for the 
numerical results. The colour indices of quarks and anti-quarks are denoted by 
and . . .} respectively. Gluons will be denoted by g with adjoint colour indices labelled 
by a, b, c, etc. 

It will be convenient to consider also fermions in the adjoint representation of the gauge 
group. These will be denoted by capital letters Q and Q and, in analogy to the fundamental 
fermions above, by Qi and Qi with an explicit flavour label. The distinction of fermions and 
anti-fermions is then induced from the flavour charge. As above the fermions Q and Q are 
understood as flavour one states, serving as a shorthand notation for Qi and Qi- 

B. Colour algebra 

The aim of the colour decomposition is to disentangle colour and kinematic degrees 
of freedom. This leads to a reorganization of scattering amplitudes in terms of ordered 
subamplitudes, which are typically easier to compute than the full amplitude itself. In 
addition, colour decomposition gives control over the interplay of colour and kinematic 
degrees of freedom. This can be used, for example, for the expansion of scattering amplitudes 
in terms of powers of the rank, Nj,, of the gauge group SU{Nc). Such an expansion can be 
exploited for efficiency gains p, |26| in numerical NLO computations. 

^ We will consider only amplitudes with strictly different quark flavours. The cases with identical quarks 
are easily obtained as linear combinations of these. 
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In a Feynman-diagram representation of a scattering amplitude, gluon interactions are 
proportional to the structure constants f"'^'^ of the gauge group, 

[T^,T^]=iV2r^'r, (2.1) 

where T'^ are the Lie algebra generators in the fundamental representation, normalised such 
that Tt{T°'T'') = 5°"^. The quark-gluon interactions carry the fundamental representation 
generators (T")/. 

We seek a uniform way of treating adjoint and fundamental colour indices, allowing the 
identification of common group theory factors. To this end we re-write all colour factors in 
terms of fundamental representation generators (T")/ using 

^abc ^ — ^ Tr([T^ T*] T'^) . (2.2) 
v2 

Contracted adjoint indices can be reduced using the Fierz identity, 

E (^")^'(^")/ = ^k^i - (2-3) 

All group theory factors can then be written in a canonical way in terms of combina- 
tions of basic group theory data: powers of N^., traces Tr (T'^^T^^ . . .T*^*^), generator strings 
(^j^airpa2 _ _ _T"-'')^ and Kronecker deltas S^. No repeated adjoint indices appear, as they are 
all reduced using the Fierz identity (12. 3p . 

Within generic formulae, when labels take boundary values, it is useful to adopt a slightly 
unconventional notation to keep expressions simple. A typical case that appears is a colour 
structure, (T"-^T"'^ . . .T"'-')/ with k ^ 0. The natural interpretation for this is to replace 
this colour structure by a Kronecker delta, 

{T'^'T'^K . .T"^)^ ^ 61 for k = 0. (2.4) 

Similarly, two boundary cases appear for traces of SU (Nc) generators, 

TriT^'T"^ . . . T"'') TriT"') = for k = l, (2.5) 
Tr(T"iT"^..r'''=) ^1 for k = 0. (2.6) 

One might find the replacement Tr(T"^T'^^ . . .T"''=) — )■ Tr(l) = Nc more natural for k = 0, 
however, we prefer to absorb factors of Nc into the definition of other objects, i.e. the partial 
amplitudes as defined below in section 111 CI 

The above replacement rules (12. 4p . (12. 5p and (12. 6 p will be understood implicitly through- 
out the present article. 



C. Partial amplitudes 



Grouping terms of scattering amplitudes according to their colour structure gives their 
decomposition into partial amplitudes. The latter are in one-to-one correspondence with 
the distinct colour structures that appear in a scattering amplitude. In this section we spell 
out the form of two-quark, four-quark and six-quark amplitudes, as they will be used in the 
later parts of the article. Wc consider these amplitudes for distinct quark flavours. The 
cases with identical quark flavours can be obtained by appropriate (anti-) symmetrization 
of quark labels and we will not state them explicitly. 

The partial amplitudes associated with the generic colour structures, 

Tr(T"^ • • • T"'^~^) (T"^+^ • • • T'S+fc-ij j^j^s+fc+a . . . j^aj+k+i^^j+k+i+i _ _ _ (2 7) 

will be denoted by, 

An;j,k,l,M, ...,j-l;jq,j + l,..., {j + ] {j + k + l)g„j + k + 2, . . . , {j + k + I + ;...), 

(2.8) 

mimicking the index structure of the colour traces and generator strings. As usual, legs 
without particle subscripts are taken to be gluons. The variables fi take values within the 
set {1, • • • , n/} according to the flavour arrangement of the quarks. At tree level fewer colour 
structures appear; the single colour trace is absent if at least one quark pair is present. The 
indices j and k,l, . . . specify the size of colour traces and generator strings, respectively. 

Partial amplitudes have symmetry properties implied by the colour structures. This 
includes a cyclic symmetry under rotation of gluons in a colour trace, 

A„.,-...(l, 2, . . . , J - 1; . . .) = ^„;,-...(2, . . . , J - 1, 1; . . .) , (2.9) 

where we rotated the gluon indices associated with the colour trace. The full symmetry 
group generated by this operation is Zj_i. Similarly, the order in which the colour structures, 
(T"™ • • •T"")j appear leaves the value of the partial amplitudes invariant, 

An-j,k,i,-(- ■■■Jq,---,{3+ k)qf^ ] {j + k + . . . , {j + k + 1 + 1)^^^ ;...)= (2.10) 

An;j,l,k,...{- ■■■,{j + k + l)g2, . . . , (j + + / + l)q^^;jq, ■ ■ ■ , (j + k)^^^] . . .) , 

and similarly for the exchange of other combinations of quark pairs. 

Explicit decomposition of scattering amplitudes into partial amplitudes will be given in 
the following. 



1. Two-quark partial amplitudes 

At tree level the two-quark QCD amplitude Q, [is], I5ll-l54j| is 



^r(l.,2^,3,...,n)= {T^^''^---T^-'-%'AT%l,,a{^),....a{n),2^), (2.11) 

where Sn-2 denotes the permutation group of n — 2 gluon labels with the amplitudes A^^^ 
the tree level partial amplitudes. The elements of the permutation group a are used in a 
two-fold way. On the one hand as a permuted list of gluon labels and, on the other hand, as 
functions, o"(A;), specifying the map of a given gluon, here gluon fc, under the permutation 
(J. The colour structures here consist of a string of fundamental generators, terminated 
with the indices of the quark and anti-quark. We suppress the coupling constants here and 
throughout this paper. 

At one-loop level the colour decomposition into partial amplitudes was given for all gluon 
multiplicities in ref. jl8|. In addition to the tree-level colour structures, colour traces appear, 

A(l„2^,3,...,n) = (2.12) 

n— 1 

Tr(T'''"(3) . . . T^-'tJ+i)) (T"''(J+2) . . . T""^'"))-^^ 

j = l (Tg5„-2/5(„_2:j) 

X (a(3), . . . , a(j + 1); 1,, a(j + 2), . . . , a{n), 2^). 

The sum in eq. f l2.12p runs over all distinct colour structures. This is achieved by a double 
sum. The outer sum runs over the group theory structures and, within each of these, the 
inner sum runs over independent gluon orderings. Given the cyclic symmetry of the colour 
traces Tr(T"'^(3) . . . T^-^tJ+i)), the full permutation group Sn~2 of the n — 2 gluons is reduced by 
the cyclic subgroups S(^n-2;j) = that leave the traces in fl2.12p invariant. We adopt here 
the notation in eqs. (12. 4p and (12. 5p for the boundary cases j = n — 1 and j = 1, respectively. 



2. Four-quark partial amplitudes 

We can follow analogous steps for amplitudes with two pairs of (distinct-flavour) quarks. 
Starting from a Feynman-diagram representation of a (gi qiq2Q29 ■ ■ ■ d) amplitude, repeated 
application of eqs. (12. 2 p and (12. 3 p leads to the colour decomposition in terms of partial 
amplitudes. 
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It is convenient to give the explicit form of the tree amphtudes for four quarks and n 



gluons 



15 



51 



54 



55|, 



(Igi) 3^2, 4^^; 5, . . . , n) — (2-13) 

n— 3 

(^T"<t(5) . . . ya^Cfc+S)^ «2 ^^J^a^^k+A) . . . J'O-aCn)^* 

k=l ■K&P2 a£Sn-4 

X A^7(7r(l,J, a{5), ...,a{k + 3), 2^^; 7r(3,J, a{k + 4), . . . , a(n), 4^^. 



For the boundary cases k = 1 and k = n — 3 the colour structures (T^'^f^) . . . 2"a„(fc+3))^*2^^^ ^^(1 
(2^a<T(fc+4) . . . T"'^("))/j'^^^ have to be replaced by 5^*^^^^ and 5^**^^ respectively (see section IXTAI) . 
P2 stands for the permutation group of the quarks {lg-^,3g^}. In an overloaded notation, 
the elements n E P2 are also interpreted as maps on the quarks and quark labels, e.g. the 
notation 7r(lqJ and 7r(l) is used to denote mappings of quark and quark label, respectively. 
The signs keeping track of the sign in the (— l/Ai'c)-term (see eq. (12. 3p ) in the colour-octet 
projector of the gluon propagator are pulled into the definition of the partial amplitudes in 
terms of colour-ordered tree amplitudes. Details can be found in the supplied data files as 
described in section IIVBI 



At one- loop level a similar colour decomposition into partial amplitudes gives 45l-l47l| . 



'^nllgn 2qj^, 3^2, 4^^; 5, . . . , ri) — 

^ ^ ^ Tr(T''''(5) . . .7^«ao+3)) 

j=l,n-3 7reP2 (Te5'„_4/S'(„_4.j) 
fc=l,n-(j+2) 

X (T"'^(J+''' ■ ■ ■ T"'-y(j+-i+k)'^^'^ (j'o.au+s+k) . . . ya<^(u)~)^«4 (2.14) 

xA„y,fc(cr(5,i+3); 7r(lgJ, (T(j+4j+2+fc), 2g^; 7r(3q2), 0'(j+3+fc,n), 452) , 

where we again use the conventions for the boundary cases given in section III A[ We also 
introduced the abbreviation (T(jj) = {a{i),a{i + 1), . . . ,a{j — l),a{j)}. 



3. Six- quark partial amplitudes 

We now consider amplitudes with three distinct quark flavours. Starting from a 
Feynman- diagram representation of a (q'l ^1 ^72^2 Q'3^3 fl' ■ • • fl') amplitude, repeated apphca- 
tion of eqs. (12. 2p and (12. 3p once again leads to the colour decomposition in terms of partial 
amplitudes. For simplicity of exposition we give the six-quark and seven-quark one-gluon 
partial amplitudes. 
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The explicit form of the tree amphtudes for six quarks is given by 15|, |5ll |5J] , 



(1,., 2^^, 3,2, 4^^, 5,3, 6-J = J2 (2.15) 
X Al;^^, (vr ( 1, J , 2,^ ; vr (3, J , 4^^ ; TT (5,3 ) , ) , 

with P3 the permutation group of the quarks {1,^,3,2,5,3}. 

For one-loop amplitudes a similar colour decomposition into partial amplitudes gives, 

3,2, 4,„ 5,3, 6,3) = J2 C(.,C(3)€(5) (2-16) 
X A-i 1 1 ( 

The generalization of the above colour decompositions f l2.16p to seven-point partial am- 
plitudes is given by, 

A(l..,2,„3,2,4,„5,3,6,3;7) = J2 H C^C^ (2-17) 

X A7;i,2,i(7r(l,J, 7, a{2gj; 7r{3,,), a{A-J; n{5,,), a{6gj) . 

The cyclic group Z3 is understood to rotate the three quarks {1,^, 3,3, 5,g} in order to 
guarantee that the gluon is emitted from all quark sources. Refined definitions of various 
elements in eq. (12.171) follow conventions in eq. (I2.15p . 

Tree amplitudes take an analogous form, as was the case for the 6-point formulas (I2.16P 
and (I2.15p . so we do not state the tree-level formula here explicitly. 

D. Primitive amplitudes 

Partial amplitudes can be further expressed as linear combinations of 'primitive am- 



plitudes' [l8|. In fact, primitive amplitudes are defined as a set of gauge- invariant, colour- 
ordered building blocks which suffice to express a given set of partial amplitudes. We discuss 
in this subsection how to characterize and generate primitive amplitudes. Our algorithm for 
expressing partial amplitudes in terms of primitive amplitudes will be given in section lllli 
Primitive amplitudes can be specified by sets of adjoint representation colour-ordered 



Feynman diagrams 18|, |49|]. These diagrams differ from the standard Feynman diagrams in 
that colour labels are dropped. At gluon vertices, for example, the structure constants /"^'^ 
are omitted. Gluon vertices dressed with colour labels are symmetric under the exchange of 
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(a) (6) 




(e) (/) 

FIG. 1: Parent diagrams with a distinct routing of quark lines: (a) a parent with a 'left-turner' 
fermion labelled by 'L', (b) a 'right-turner' fermion, labelled by 'i?', (c) 'n^-left-turner' (rifL) 
fermions, (d) 'n/-right-turner' {ufR) fermions. Two routing labels are associated to the distinct 
quark flavours Q and Q2'- (e) 'left-left-turner' {LL) fermions, (f) 'left-right-turner' {LR) fermions. 
Note that in (f) the first fermion does not in fact enter the loop. In general as many L/i?-labels 
appear as there are distinct external quark flavours. Quarks are associated with capital letter 
Q's, as in primitive amplitudes they arise from adjoint representation colour assignments in our 
algorithm. 
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legs. Due to the anti-symmetry of the structure constants f , the remaining kinematic part 
is necessarily anti-symmetric under exchange of legs. In colour-ordered Feynman diagrams 
one must preserve the ordering of the external legs to keep track of the signs associated 
with the vertices. Given a set of colour-ordered Feynman diagrams one uses colour-ordered 
Feynman rules (see e.g. |48;]) to obtain amplitudes and eventually primitive amplitudes. 

We require the set of all primitive amplitudes needed to compute a given scattering 
amplitude. In particular, we will need their representation in terms of sets of Feynman 
diagrams. This may be generated systematically in two steps, which we will describe in full 
below: 

1. generate 'coarse' primitive amplitudes, 

2. split the coarse primitives according to the routing of fermions around the loop. 

The coarse primitives can be thought of as sums of planar, ordered diagrams, obtained using 
colour-ordered Feynman rules. Each coarse primitive amplitude corresponds to a particular 
ordering of the external legs. 

In practise, we obtain the coarse primitive amplitudes starting from a colour-dressed 
representation of the amplitude by taking all partons, including quarks, in the adjoint rep- 
resentation of the gauge group. Performing all the colour algebra reductions, one finds a set 
of terms, some with single and some with double traces of fundamental generators, just as 
is the case for purely gluonic amplitudes. The coarse primitive amplitudes 



Secondly, a finer splitup of the coarse primitive amplitudes into actual primitive ampli- 
tudes is obtained by analysing the fermions' routing around the loop. This can be considered 
a generalisation of the familiar case of diagrams with a fermion circulating in the loop not 
mixing under gauge transformations with diagrams with a gluon in the loop. Similarly, 
fermion routing information in Feynman diagrams allows to split coarse primitive ampli- 
tudes into finer gauge-invariant subsets. The routing information of a fermion line specifies 
whether the fermion turns left or right upon entering the loop. In the case where the fermion 





are then identified as the coefficients of the single trace structures 




(2.19) 
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does not enter the loop, the left /right label is assigned according to which side of the loop 
the fermion passes. We follow the fermion line starting from the anti-quark. In this way, 
each fermion line is labelled as a left- or right-turner, denoted by capital labels L and R, 
respectively. Some examples can be found in fig. [T] 

Algorithmically, the routing information of a colour-ordered Feynman diagram can be 
determined by cutting a single loop propagator. If no gluon propagator is available, the 
Feynman diagram has to correspond to an term. We then cut the internal fermion line. If 
a gluon loop-propagator is available we instead cut that one. After cutting the propagator the 
Feynman diagram corresponds to a diagram of a colour-ordered tree amplitude. The fermion 
routing information is then in one-to-one correspondence to the ordering of external fermions 
relative to their anti-fermions and may be read off directly. The coarse primitive amplitudes 
can then be sorted diagram by diagram into the finer classes of primitive amplitudes with 
definite fermion flow. 

The explanation of the gauge invariance of the primitive amplitudes as defined above has 
been given in ref. 18|. The basic reasoning relies on a non-standard colour charge assignment 
to quarks and gluons. Primitive amplitudes of standard QCD can typically be interpreted as 
partial amplitudes of a theory with non-standard gauge groups and colour charges. Gauge 
invariance of a class of Feynman amplitudes is then ensured from the associated gauge invari- 
ance of partial amplitudes. The non-standard colour assignments include various bifunda- 
mental and adjoint representations in product gauge groups SU{Ni) x SU{N2) x SU^Ns) . . . 
for quarks and gluons respectively. The basic logic is only a minor generalisation of that 
presented in the original literature {isl], to which we refer for further details. 



1. Notation 



Pictorially, primitive amplitudes can be characterized by a small set of 'parent' Feynman 
diagrams. These are the representative colour-ordered Feynman diagrams of a primitive 
amplitude with the maximal number of loop propagators. An additional requirement is 
that the complete set of Feynman diagrams associated with a primitive amplitude can be 
obtained from the parent diagrams through pinching of propagators and 'pulling out' of tree 
amplitudes. 

Example parent diagrams of six distinct primitive amplitudes are shown in fig. [T] The 
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diagrams are representative of classes of colour-ordered Feynman diagrams. Parent diagrams 
(a), (b), (c) and (d) have identical external states, as do diagrams (e) and (f). They differ 
by the routing of the fermion lines around the loop and are gauge invariant individually^. 
The routing information of the fermions is given in terms of left- and right-turner labels; L 
and R. Closed fermion loops are specified by an additional label nf. Such routing labels 
must be specified for each fermion line. 

So instead of specifying a primitive amplitude by the set of colour-ordered Feynman 
diagrams contributing to it, we can be more concise, and specify only 1) an ordered set of 
external states and 2) routing data of fermion lines. 

We will use this information to specify primitive amplitudes. In our conventions, for the 
amplitudes shown in fig. [T] we will use the following notation. 



(a): 


A^iQ,---g,9,Q,g,9,9,g,---) 


(6): 


A^iQ, ■ --g^g^Q^g^g^g^g, ■ ■ ■) 


(c): 


A''f^{Q,...g,g,Q,g,g,g,g,.. 


(rf): 


A''f^{Q,...g,g,Q,g,g,g,g,.. 


(e): 


A^^{Q,Q2,g,Q2,- ■ ■g,g,Q,g, 


(/): 


A^^{Q,Q2,g,Q2,---g,g,Q,g, 



Fermion lines will be given ascending flavour numbers starting from 1 (we identify Q = Qi). 
The routing information will be specified by a list of labels L, R. L stands for a fermion 
that turns left upon approaching the loop, while R stands for a right turning fermion. The 
first label in the list corresponds to the first flavour-line, the second label to the second line 
etc. Closed fermion loops are indicated by an additional label n/. External states are given 
in clockwise order and we exploit the rotational symmetry (described in the next section) 
to rotate quarks with flavour one into the flrst position. 

2. Symmetry properties 

Pure-QCD primitive amplitudes enjoy a variety of symmetries. These are useful for 
reducing the number of independent primitive amplitudes that need to be computed. We 

^ More precisely, the primitive amplitudes they correspond to are gauge-invariant. 
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use the symmetries below to write primitives in a standard form starting with a left-turner 
label, with the quark of flavour one in the first position of the particle labels, 

A^-(Q,...) or A"/^-(Q,...). (2.20) 

The relevant symmetry transformations are described below. In addition, an extended set 
of relations between primitive amplitudes will be discussed below in 13.101 

First of all, cyclic rotation of the labels does not change the primitive amplitudes, 

A-(l,2,...,n) = A-(2,...,n,l). (2.21) 

Another symmetry arises from flipping over the colour-ordered diagrams of a given primi- 
tive. This reverses the ordering of the external states and exchanges all left-and right-turners 
up to a sign. In the case of a two quark line amplitude we have, 

A^^(l„ 2, fc^, mg^, n) = (2.22) 

(—1)"" A^^{lq, ra, . . . , rn-q^, . . . ,lq2, . . . , kq, . . . , 2) . 

A further generic relation between primitive amplitudes originates in the reversal of a 
fermion's arrow. This leads to a flip of the routing label, 

^""^'■■■(. . . , kq^, . . . , /q^, . . . , n) = ^^■■■■^'•••(^. . . , /cq., . . . , /q., . . . , n). (2.23) 

This property is naturally generalised to the case with multiple quark lines. Each quark 
line may be reversed and the associated turner label flipped, L ^ R. We will not use this 
transformation here, preferring to keep the distinction between quarks and anti-quarks. 



E. Colour decomposition of partial amplitudes. 

In this section we present our notation for partial amplitudes given as linear combinations 
of primitives. As the expressions are lengthy we prefer to write them in a simplified form, 
which is also more convenient for use with a computer program. We display only a small 
sample in the main body of this paper, and attach the remainder as text files j56| . Additional 
information concerning the supplied data-files can be found below in section lIVBi 

To explain our notation we refer to fig. [2] where we show our result for the four quark 
amplitude. The line starting with '* Partial' indicates the beginning of the expression 
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and is followed by a specification of which particular partial amplitude is being considered. 
The notation for this mimics the argument of the partial amplitudes (12. 8p . It takes the 
form of a comma separated list with each element of the list denoting an SU{Nc) colour 
structure. Such a colour structure is unambiguously specified by a list of particles. Quarks 
in the fundamental representation are denoted by q, and may have an additional numerical 
label to distinguish different flavours (g, q2, q3 ...). Anti-quarks are written qb, qb2, qb3 
etc. Momentum labels appear in parentheses, e.g. denotes a quark that carries the 

momentum ki. For the case of fig. |2], the colour structure is Sl'^Sl'^- 

There follows specification of the born contribution to this partial amplitude in terms of 
colour-ordered tree-level amplitudes. The line starting with the label, 

born {. . . I . . . 

etc. 
} 

shows a list of particles in the adjoint representation and, separated by a vertical line, '|', the 
coefficient of the given tree within the partial amplitude. Labels with capital letters Q, Q2 
etc. are used for these kind of quarks. For partial amplitudes that are a linear combination 
of colour-ordered tree amplitudes, multiple line entries are used, each displaying a list of 
particle labels and the weight of the associated tree amplitude within the partial amplitude. 
In more conventional form the four quark born amplitude is given by, 

ATH,, 2„ 4-J = A'-%1q, 2q, 3g„ 4qJ x (^6]^ - i- S^^S]^ . (2.24) 

* Partial: { q(l) qb(2) , q2(3) qb2{4) } 

born { Q(l) Qb(2) Q2(3) Qb2(4) | (-1/Nc) } 
loop { 

LL I Q(l) Qb2(4) Q2(3) Qb ( 2 ) | 1/Nc"2 
LR I Q{1) Q2{3) Qb2(4) Qb ( 2 ) | l/Nc'^2 + 1 
LR I Q(l) Qb2(4) Q2(3) Qb ( 2 ) | 1/Nc"2 
LL I Q{1) Qb{2) Qb2(4) Q2(3) | -1/Nc^2 
nfLL I Q(l) Qb2(4) Q2(3) Qb ( 2 ) | -nf/Nc 

} 

FIG. 2: Four quark partial amplitude decomposed into primitive amplitudes, in the notation of 
the attached file. 
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The displayed partial amplitude corresponds to the second term in eq. fl2.24p and the ordered 
born amplitude is A*'^^'^{1q, 2q, Sqj, 4q^) weighted by a factor of — 

The remaining entries specify the one-loop partial amplitude in terms of primitives. This 
entry starts with the header '/oop', 

loop { 



etc. 

} 

with the individual lines each specifying the turner labels, an ordered list of external states 
and the weight of the primitive amplitude. In total this represents a sum of primitive 
amplitudes, each with a Nc- and/or nj-dependent coefficient. For the present example, 
displayed in fig. [21 the partial amplitude is given by, 

A(l„25,3,2,4gJ = ...+ /l4;i,i(l„25;3,2,4gJ5:j5g + ... , (2.25) 

+ W ^^""^IQ' 2q) - 4 3qJ 

-'■^ c -^^c 

-^^"^^^(1q,4^^,3q2,2^). (2.26) 

The five primitive amplitudes appear in the same order as in fig. |2l where they are given 
on separate lines. Each primitive amplitude is specified by a list of particles in the adjoint 
representation and string of fermion left /right-turner labels. We denote adjoint quarks with 
a capital Q. Each quark line has a 'turner label' L or i?, which described the direction it 
turns upon approaching the loop, as described in section III D[ 

The above expression can be compared to the one given for the partial amplitude denoted 
by in eq.(2.15) of ref. {iq!]. To simplify this comparison we set the number of scalars 
rZg = and drop the top-quark contributions in the expression for Aq.2. We can omit the 
colour-neutral lepton states, as they play no role in this colour decomposition. To match 
the expression for Aq.2 with our eq. fl2.26p we have to replace 2 -H- 4 in the expression for 
Aq,2- We can then identify ^+"(1,3,4,2) = A^^(l, 3, 4, 2), A++(l,4,3,2) = A^^(l, 4, 3, 2), 
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A^'(4,3,l,2) = A^^(l,2,4,3) - A^^(l,4,3,2) and 4, 3, 2) + 4, 3, 2) = 



— A"^^^(l, 4, 3, 2). (Compared to [l9| we differ here by an overall sign in the definition 



of the nj-terms, which introduces this apparent relative minus-sign.) This identification can 
also be inferred from comparison of parent diagrams. The relative minus sign in the expres- 
sion for A^''{A, 3, 1, 2) appears from the asymmetry of the fermion-fermion-gluon three-point 
vertex. Upon reversing the order of the quarks Iq and 2q and rotation of the labels one finds 
the identity, -A^^{1, 2, 4, 3) = A^^{2, 1, 4, 3) = A^^{1, 4, 3, 2), with the relative minus sign 
absorbed into the re-ordering of external legs. 

With this comparison we conclude the explanation of our conventions. 

III. COLOUR DECOMPOSITION ALGORITHM 

In this section we will describe our setup for obtaining partial amplitudes in terms of 
primitive amplitudes, for arbitrary processes. The algorithm we use is based on analysing 
Feynman diagrams using their colour information. We identify partial amplitudes in terms 
of linear combinations of Feynman diagrams. Similarly, we express primitive amplitudes in 
terms of Feynman diagrams. Finally, we express partial amplitudes in terms of primitive 
amplitudes by solving the linear set of equations to eliminate the explicit dependence on 
the diagrams. In addition to the colour decomposition of the partial amplitudes, we find 
non-trivial relations between primitive amplitudes, which arise due to the redundancy of the 
linear equations we solve. 



We begin by generating all one-loop Feynman diagrams using QGRAF [57|]. We remove 
diagrams with contact terms (four-point gluon vertices), tadpoles and those with bubbles 
on external lines. Let us suppose there are A^^^ diagrams after this procedure. ^he output of 
this program is then processed using the computer algebra package FORM [58]. Using this 
package we dress the diagrams with colour - gluons in the adjoint and quarks in the funda- 
mental representation of the gauge group - and simplify the colour algebra using eqs. (12. 2p 
and (12. 3p . There is no need to substitute the kinematic parts of the Feynman rules, so we 
leave this information implicit. This procedure gives us an expression for the amplitude in 




A. Setup 
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terms of sums of diagrams multiplied by the various colour structures that may appear, 

Np 

^(K,iJ) = X]Ar"^^C.({a„iJ). (3.1) 

i=l 

On the left hand side is the fully colour-dressed amplitude, which depends on the set of 
colour charges {ag, iq} of the external particles. The right hand side is its decomposition in 
terms of Np partial amplitudes ^p'^''*"^^^ multiplying colour structures Cj. Partial amplitudes 
and their colour structures are discussed further in section III CI Here and in the following 
we leave implicit all dependence on momenta and helicities. 

After processing the diagrams in this way, we find expressions for the ^p^''*"^^ as sums of 
diagrams, 

^partial for l = l,...,Np. (3.2) 

i=i 

Here Ki^ is a Np x N^ matrix which describes how the j'th diagram contributes to the i'th 
partial amplitude. The entries of the matrix Ki ^ consist of simple functions of the rank of 
the gauge group, N^. The di are the kinematic parts of the set of Feynman diagrams, but 
we will not need their explicit form here. It is sufficient to think of them as tags for the 
individual Feynman diagrams. 

With this in hand, we now discuss how the primitive amplitudes are expressed in terms 
of diagrams. We first go back and dress the original one-loop Feynman diagrams with colour 
assuming that all external particles, including quarks, live in the adjoint representation of 
the gauge group. Upon simplifying the colour algebra, we find two types of term — those 
with single traces, and those with double traces, exactly as one finds for purely gluonic 
amplitudes, 

_^adjoint(^|^^^ aj) = J2 Tr(T''-(i)T''-(2) ■ ■ ■ T'^-(")) A''°''''°{a) + double trace terms. 

(3.3) 

In this equation we have added the label 'adjoint' to remind the reader that this is the 
amplitude with all external particles in the adjoint representation. The sum is over the 
permutation group S'„ of n elements modulo cyclic rotations Z„. Some permutations may 
not appear when fermion lines of the diagrams cannot be arranged in a planar way. We 
drop all double-trace terms. Note that there is no approximation being made here — this 
is simply the definition of the primitive amplitudes. (The double trace contributions can be 
fully reconstructed from the primitive amplitudes, though we will not need to do this here.) 
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If we were to restrict ourselves to the case of amplitudes with only adjoint particles, then 
we could now identify the coefficients A'^"^'^^'^{1,2, . . . ,n) of the single trace structures as 
primitive amplitudes, 

^---(afc) = k = l,..., iV™ , (3.4) 

i=i 

where {L^"^'^^'^)k^ is a coefficient matrix with integer entries. 



However, as described in section IIIDl following ref. 18| for the case of a single fermion 
line, when fundamental quarks are present these objects can be split further into finer gauge 
invariant pieces. These pieces are distinguished by the direction of the fermion line relative to 
the loop. Following the fermion line into the diagram, one defines 'left-turner' (L) fermions 
as those which pass to the left of the loop, while 'right-turner' (R) fermions pass to the 
right. If the loop is in fact a closed fermion loop, the diagram acquires an additional label 
Uf. These concepts have already been discussed in section HTDl We extend this grouping of 
diagrams to generic processes by labelling each fermion line as a left- or right- turner. We 
thus take the primitive amplitudes to be 

A^i-^^«(l,...,n), (3.5) 

where Nq is the number of fermion lines and Di G {L, R} labels each as either a left or 
right turner. In the argument list we have suppressed the identities of the 2Nq fermions and 
n — 2Nq gluons. 

Let us suppose there are A^prim primitive amplitudes. Using our setup we derive an 
expression for each of them as a sum of diagrams, just as we did for the partial amplitudes 
in eq. I K2^ . 

Nd 

Af' = (ak) = Y,Lk'dj, k = l,..., iVpHm (3.6) 
i=i 

where the direction labels are implied. Here Lk ^ is a A^prim x matrix describing how the 
j'th diagram contributes to the k'th primitive. 

We now attempt to express the partial amplitudes of eq. fl3.2p as linear combinations of 
the primitive amplitudes of eq. ( 13. 6p . We write 

prim 

^partial ^ ^ (3.7) 
fc=l 
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where i runs from 1 to Np, and look for solutions for the coefficients Zi ^ of this set of linear 
equations, 



Zi ^ Lu ^ 



(3.8) 



The solutions of this set of equations are typically not unique as will be discussed further in 
section IIIIBi 

We have confirmed that the solutions match those obtained from the known all- 



multiplicity expressions in ref. [17| in the case of purely gluonic processes, and ref. [18 1 
for processes with one quark line. For processes with more than one quark line we repro 
duce the closed form analytic expressions available in refs. 
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46| and [47[. Beyond that 



we are able to present new results as presented in more detail in section IIV[ 



B. Relations 



The linear equations (13. 8p typically give rise to non-trivial solutions {{Z 



(3.9) 



This means that the Feynman diagrams associated with a linear combination of primitive 
amplitudes add up to zero. Interpreted for the primitive amplitudes this implies linear 
relations, 

^ "prim 

= ^ (Zf)'' ^ = 1, • • • , AT-i'^*'"- , (3.10) 



fe=i 

for the case that A^''^' solutions {Z^'^^)^ have been found. In order to verify the relations no 
detailed knowledge of the Feynman amplitudes is necessary. They arise at the diagrammatic 
level when taking into account the anti-symmetry of colour-ordered three-point vertices. 
It is tempting to relate the above relations to the ones that are known for colour-ordered 



tree amplitudes, such as the U{1) decoupling identity 
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52 



izations 



21| through the Kleiss-Kuijf relations 59j. The role of these re 



53| and non-abelian general- 



colour-decomposition has already been pointed out some time ago 
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ations for loop-level 
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60| and recently 



reviewed in ref. SOj. We will not attempt a detailed analysis into this direction but make 



some general remarks about the observed relations (I3.10p . 

One remark concerns the importance of having multiple fermion lines to obtain the above 
relations. When considering fermion lines the anti-symmetry of the fermion-fermion-gluon 
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three point vertex implies certain fermionic amplitudes may be related. The complete list 
of all relations is a by-product of our way to obtain the colour decompositions of partial 
amplitudes. By inspection we observed that the relations in eq. (13.1 Op originate from this 
anti-symmetry. We will discuss two examples below in section IIII B 1[ As the multiplicity of 
quark amplitudes is increased by addition of gluons, the number of relations increases due 
to the different ways gluons can be added to ordered amplitudes. The new relations appear 
as descendants of the ones of the purely fermionic amplitudes. 

A second remark is about an application of the above relations. For numerical evaluation 
of scattering amplitudes (see later in section lIVp we exploit the relations between primitive 
amplitudes (I3.10p to optimise caching and thus run-times. In the explicit expressions given, 
we pick methodically a subset of primitive amplitudes using the relations. In fact we sort 
primitive amplitudes according to the number of propagators of their parent diagrams. We 
then apply the relations to eliminate primitive amplitudes starting from the ones with the 
minimal number of parent propagators as part of the loop. This procedure minimizes the 
number of primitive amplitudes that need to be computed. However, it does not take into 
account which primitive amplitudes are most easily computed. 

Similarly, one could make use mostly of primitive amplitudes with a minimal number of 
parent loop propagators. For a unitarity-based algorithm this choice would be beneficial, 
reducing the number of unitarity-cuts to consider. We leave a more thorough analysis of the 
relations between primitive amplitudes to the future. It seems likely that an understanding 
of the relations between primitive amplitudes (IS.lOp will be helpful towards establishing all-n 
formulae for the colour decompositions, which are beyond the scope of the present article. 

1. Example Relations 

For the set of primitive amplitudes we are dealing with the origin of these relations are the 
symmetry properties of the quark-quark-gluon vertices. A simple example with an internal 
fermion line is given by, 

A°^^L(1q, 5q^, 4q„ 2, 3q)+A"^l^(1q, 4q„ 5q^, 2, + A^^'^'^ilQ, 5^^, 2, 4q„ 3q) = , (3.11) 
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FIG. 3: Parent diagram representation of relations between multi-quark primitive amplitudes: 
Diagrams (a) show a vanishing sum of primitive amplitudes associated to parent diagrams with a 
closed fermion loop, as in eq. (13. lip . Diagrams (b) show a vanishing sum of primitive amplitudes 
with at least one gluon in the loop, the pictorial form of eq. (|3.12p . The relations originate from 
exchanging the ordering of the fermion pair {4:q,^ , 5q^}, which gives a relative sign. Parent diagrams 
represent classes of colour-ordered Feynman diagrams. The contributions of all colour-ordered 
Feynman diagrams cancel once the asymmetry of the fermion-fermion-gluon vertices is taken into 
account. 

with the parent diagrams shown in fig. [31 Similarly, relations my be found for diagrams with 
no closed internal fermion line, 

ALL(1q, 3q, 5q^, 4q„ 2) + 3q, 4q„ 5q^, 2) + ^^^(Iq, 3q, 5q^, 2, 4qJ = , (3.12) 

with the parent diagrams shown in fig. [3] (b). 

For both diagrammatic equations, (a) and (b) in fig. |3l the same mechanism is at work. 
Contributions from the first two diagrams, respectively, with a direct (4Q2,5Q_^,gluon)-vertex 
cancel against the ones with a direct (5Q^,4Q2,gluon)-vertex, given that the exchange of 
quarks Aq^ ^ introduces a relative minus-sign. The remaining contributions include 
Feynman diagrams with the gluon moved onto the fermion line. For the first parent diagram 
in in fig.|3]that would be gluon emission from 4q^, and for the second parent diagram emission 
from 5q^. These contributions are distinct and it is the role of the Feynman diagrams 
represented by the third parent diagram, to cancel off these remaining pieces. 
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FIG. 4: Feynman diagrams di contributing to the four-quark amplitude. For simplicity we consider 
only box and triangle diagrams. These are sufficient to obtain the colour decomposition of the four- 
quark partial amplitudes. 

C. A four-point example 

To illustrate the method described in the previous section we present here a fully worked 
example. We consider the amplitude for scattering of two pairs of different flavour quarks, 

A{lg,2-,,3,„A-J- (3.13) 

To keep things as simple as possible we will ignore contributions with closed fermion 
loops, though in our final results these pieces are of course included. We make one further 
simplification for transparency of the example: we consider six Feynman diagrams (A^^^ = 
6), illustrated in figure IH and drop bubble diagrams. These diagrams are sufficient for 
the present example. After dressing them with colour (quarks and anti-quarks live in the 
fundamental and anti-fundamental representations of SU{Nc) respectively), and simplifying 
using eq. (EH]), we find two distinct colour structures {Np = 2), 

Ci = 51J<5g, C2 = 61^61. (3.14) 
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The amplitude is a linear combination of Ci and C2, 

-^(Ig; 2g, 3^2, 45^) = A\^'^^^^^{lq, 2g; 3,2, 4qJ Ci + Al^'^^^'^^{lg, 4^^; 3^2, 2g) C2 ■ 



(3.15) 



The expression of the partial amplitudes, A^ (Ig, 2^; 3^2? 4^^) and A2 (1^, 4^^ j 3^2, 2g), 
in terms of the Feynman diagrams in figure H] is described in (13. 2p through the 2x6 matrix 
of coefficients Ki\ 



-1 



1 

Af2 



-1 



1 

Ar2 



1 

Af2 



2 



1 

iV2 



AT - A. 



(3.16) 



when contracted with the vector of Feynman diagrams dj. The top row of this expression 
relates to Ci, and the bottom row to C2. 

To find the primitives, we must dress the amplitude with colour assuming all particles 
are in the adjoint representation, and then extract the primitives as the coefficients of the 
single trace structures. We find 







'2q'4q^'3q2) 


— —dii 


(3.17) 


A^^( 


1q 


) 3q2 5 ^Q^) 2q) 


= -d2, 


(3.18) 


j^LL 




>4q2'3q2,2q) 


= di + d3 + de, 


(3.19) 




1q 


' 2q, 3q2, 4qJ 


= di, 


(3.20) 




1q 


5 3q2 5 4q^5 2q) 


= d^ — d^ — di, 


(3.21) 


j^LR^ 


1q 


' 4q2' '^^2 5 2q) 


= d2. 


(3.22) 



Already at this level we observe that the primitive amplitudes are related as. 



A^-{1q,2oAo.,^c 



1q'3q2,4q ,2q) 



-A^^(1q,2q,3q2,4qJ 
-A^«(1q,4q^,3q2,2q) 



(3.23) 



These relations are an example of the ones discussed in section IIII B[ They are not direct 
consequences of the symmetry properties of the primitive amplitudes f l2.2ip . f l2.22p and 
f l2.23p . Rather they are explained by the anti-symmetry of the colour-ordered fermion- 
fermion-gluon vertex. In eqs. f l3.23p the exchange 3^4 and LL o LR relates the primitives 
up to a sign. This can also be inferred from the diagrams d2 and ^3 in fig. |H Following 
our earlier notation ( 13. 6 p the equations ( I3.17p -( l3^22|) are expressed through a matrix of 
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coefficients Lk^ and a vector of primitives A^'', 

/o 0-1 oo\ 
0-10 
1 1 1 
1 
-10-1010 
\ 1 / 



( A^^(1q,2q,4q^,3q, 



^'^'^(1q,3q2,4q_^,2q 
VA^«(1q,4q^,3q2,2q 



(3.24) 



We note that L is generally not square, though in this case it is because we chose a re- 
duced number of Feynman diagrams, which happens to be equal to the number of primitive 
amplitudes. 

The solution to equation (13.71) can be found by eliminating the di from eqns. (13.21) using 



le coefficient matrix Z,- 



61| 



linear combinations of the primitive amplitudes (13. 6p . We obtain t 
by solving the linear system of equations (13. 8p with Mathematica 

It turns out that there are many different solutions, reflecting the relations (13.231) between 
primitive amplitudes. This is a new feature for multiple quark line processes that is not 
observed in the purely gluonic, or single quark line cases. In fact, the matrix {L^^} is a 
6x6 matrix with rank four. The associated null space (13.91) is two-dimensional and given 
by the two vectors {Zl^^Y with i = 1, 2. In matrix notation we have. 



rc\\ k 



10 10 
1 1 



(3.25) 



The null vectors imply the same relations given already in eqs. (I3.23p . 



The solutions are represented clS cl Np X N, 



prim 



2x6 matrix Z as follows 



1 



a 1 + 



/3 
6 



(3.26) 



where a, /3, 7 and 6 are arbitrary complex numbers parameterising the solution space. The 
ffist row of Z gives the six coefficients of the primitive amplitudes, in the basis (I3.17p . needed 
to construct the coefficient of the colour structure Ci. The second row does the same thing 
for the structure C2. 

Anticipating the explicit results in section |IV] we relate the solutions (13.260 to the at- 
tached data-ffies. Compared to the explicit colour decompositions given in the attached ffie 



26 



(partials_4q.dat) we have the following choices for the parameters, a = 0, (3 = 7 = 0, 

c 

and 6 = In this way we avoid the use of two out of six primitive amplitudes. 



IV. RESULTS 



General expressions for the decomposition of partial amplitudes in terms of sets o 
itive amplitudes are available in the literature for two classes of process, the n-gluon 
and the two-quark n-gluon 18| QCD amp 



m 



)rim- 



four-quark amplitudes was given in 
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itudes. A fixed-multiplicity decomposition of 
46j and for an additional gluon in 47|. In fact 
these expression were first given including particles neutral under the gauge group. Here we 
generate sets of partial amplitudes algorithmically for a given scattering process with em- 
phasis on the cases of six- and seven-parton QCD amplitudes including four quarks and six 
quarks. The generalisation of these to amplitudes including a single colourless vector-boson 
is straightforward, and will be discussed in section IIV CI In the recent results UJ, for 
V^-boson plus four-jet NLO cross sections, a leading-colour approximation was used in the 



virtual contribution. In the following we will present a corresponding full-colour distribution. 



A. New colour decompositions for QCD amplitudes 

We present first our results for partial amplitudes in terms of primitive amplitudes with 
six and seven coloured partons. We focus entirely on processes with distinct quark flavours. 
This incurs no loss of generality, because all other cases can be derived as simple linear 
combinations of these. 

As an example we show in fig. [S] one of the partial amplitudes relevant for six-quark 
processes. In this case there are 26 contributing primitives, each with its own A^c-dependent 
coefficient. There are three turner labels, one for each quark line. As described in section UTTl 
when there are multiple quark lines present the expressions are not unique, as the primitive 
amplitudes satisfy a set of relations. 

It is interesting to consider in comparison the leading-colour part of the above partial 
amplitude, as shown in fig. O Formally, the leading colour part of the partial amplitude is 
here defined as the leading terms in the limit Ac ^ 1 and Uf ^ 1 with the ratio Nc/uf fixed. 
The number of contributing primitive amplitudes is strongly reduced when subleading in 
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Partial: { q(l) qb2(2) , q2(3) qb3(4) , q3(5) qb(6) } 
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} 

FIG. 5: A six-quark partial amplitude given as a decomposition into primitive amplitudes. The 
shown partial amplitude is an extract from the file partials_6 q.dat in the ancillary material 



* Partial: { q(l) qb2(2) , q2(3) qb3(4) , q3(5) qb(6) } 

born { 0(1) Qb2(2) Q2(3) Qb3(4) Q3(5) Qb(6) | 1 } 
loop { 

LLL I Q(l) Qb2(2) Q2(3) Qb3 ( 4 ) Q3(5) Qb(6) | Nc 
nfLLL I Q(l) Qb2(2) Q2(3) Qb3(4) Q3(5) Qb ( 6 ) | nf 



FIG. 6: Leading colour six-quark partial amplitude given as a decomposition into primitive ampli- 
tudes. 

Nc contributions are removed. This is one of the reasons for using this approximation when 
constructing NLO partonic Monte Carlo programs — the leading-colour amplitude can be 
numerically evaluated much faster. 

The leading-colour contribution to the decomposition of quark amplitudes is easily ob- 
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ia) ib) 

FIG. 7: Generic leading-colour parent diagrams for the leading partial amplitudes: (a) a {LL ■ ■ ■ L)- 
turner and (b) a {njLL ■ ■ ■ L)-turner. External gluons may be added if directly emitted from the 
gluonic loop propagators in (a) or from the fermion loop in (b). 

tained for arbitrary multiplicities following a few simple rules: 1) partial amplitudes with 
vanishing or subleading (in 1/Nc) tree- level contributions are set to zero. 2) the particle 
orderings of the leading-colour primitive amplitudes follow that of the born amplitude. 3) 
Two particular types of left-right-turner labels are singled out leading to the two types of 
parent diagrams as shown in fig. [71 In (a) we show the primitive amplitude with multiple 
left-turner quark lines including gluon propagators in the loop. It is denoted by the label 
[LLL . . .). The associated parent diagram has both gluon and fermion propagators. The 
number of propagators matches the multiplicity of the amplitude, here meant to be n. The 
second primitive amplitude, (b), has a closed fermion loop with each fermion line coupled to 
it via gluons. The label specifying the fermion routing is then {ufLLL . . .). The associated 
parent diagram has only fermion propagators. The general structure is manifest in the two-, 
four- and six-quark cases discussed here. 



B. Explicit results in attached text files 

The expressions for the partial amplitudes can be found in the location 
spective text files 

partials_nggnqq.dat (4.1) 



56| in the re- 
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provided for all independent partial amplitudes with Ug gluons and Uq quarks. The format 
of the files has been described in section III E[ We give the cases, 



{K,^.)} = {(0,4),(1,4),(2,4),(3,4),(0,6),(1,6)}. 



(4.2) 



17|, 



18| . this completes the 



Together with the (n, 0) and {n — 2, 2) assemblies given in refs. 
colour decomposition of QCD amplitudes up to and including seven colour- charged external 
states. 

We confirmed that two independent implementations of the algorithm described in this 
paper produce identical numerical results for partial amplitudes. The explicit expressions 
from each implementation match up to application of the relations described in section UlI B[ 
Furthermore, the expressions for the partial amplitudes are checked against the existing 
results in the literature where available. In particular, for all new six-parton and seven- 
parton expressions we confirmed numerically that the virtual squared matrix elements give 



the expected pole structure [62| for single and double poles in the dimensional regularization 



parameter e. We used BlackHat 



3l| and SHERPA 



for this numerical check. A 
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further numerical check is provided by the recent public release of Helac-ILoop 
which we have used to check all our 6-parton results. 

In addition to the expressions for pure QCD amplitudes, the attached files (for the down- 
load location see [56|) 

partials_nggnqq21.dat (4.3) 

give their generalisation to include a colour neutral lepton pair. These contributions are the 
ones relevant for the computation of FT + 4-jet production, as we discuss in the next section. 



C. Application to W + 4, jets 

The results we have presented are crucial for the application of the colour-ordered ap- 
proach to vector-boson production in association with multiple jets at NLO. With the colour 
decomposition of the pure QCD amplitudes one can obtain rather simply their generalisa- 
tions to include colourless states. We discuss here how to obtain the colour decomposition 
of iy-|-jets amplitudes, and give a comparison of leading-colour versus full-colour differential 
distributions of the virtual part oiW + 4-jet production at the LHC. 
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1. Inserting colour-less states. 



Amplitudes with additional colourless objects - leptons, vector bosons, Higgs bosons, etc. 
- can be accommodated in the decomposition of partial amplitudes into primitive amplitudes 
starting from the pure QCD decomposition. 

For W + n-jet production we consider the decay of the W^-boson into a lepton pair, ui 
and ^. These amplitudes can be obtained by first computing amplitudes with a virtual 
photon that is emitted from a quark q and decays to a charged lepton pair (^l,^/?). In a 
second step, the conversion to a lepton and a neutrino {vilAr) is done by multiplicative 
factors including cou plings and VT-boson propagator terms. Details about this conversion 



are explained in refs. 
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65| . With this approach we consider amplitudes with a lepton pair 



coupled to a quark pair via an off-shell photon. In fact, we do not need to consider all such 
amplitudes; we can drop the emission of the photon from an internal quark loop as these 
kind of emissions are forbidden for ly-production due to flavour non-conservation. 

With the basic understanding that the lepton pair is coupled to the quark with flavour 
one, it is straightforward to transcribe the pure QCD partial amplitude into the one including 
a virtual photon decayed into a lepton pair, as displayed in fig. [HI The full scattering matrix 
elements can be built from these basic amplitudes [l9|. 



2. Full colour differential distributions. 

The transverse momentum, p^, of the accompanying jets is one of the most important 
observables in W-hoson production. For W + 4-jet production at NLO the differential pr 



distributions of the leading four jets have been given in [ll|]. In that paper a leading-colour 
approximation was used for the virtual parts^. Here we compute for the first time the 
subleading-colour corrections and compare the size of leading-colour and full-colour virtual 
contributions. We use the same basic setup and cuts as in pjj] . The virtual contributions are 
computed using on-shell methods via the BlackHat package, while SHERPA is used to 
perform the phase-space integration. The combination of these two programs has been used 



extensively for studying hadron collider phenomenology at NLO 



11 



12 
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66 



67|, 



and is very well tested. For simplicity we will here focus on the px distribution of the fourth 



The other parts - the born and real emission - were evaluated without any colour-based approximation. 
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Partial: { q{l) qb2(2) , q2(3) qb3(4) , q3(5) qb(6) ; 1(7) lb (8) } 
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FIG. 8: Six-quark two-lepton partial amplitude given as a decomposition into primitive amplitudes. 
The leptons couple via virtual photon exchange to the quarks with flavour one, Q and Q. These 
partial amplitudes are building blocks of the scattering amplitudes for ly-boson production. 



jet. 



ll| by terms 



The leading-colour approximation used here differs from the one in ref. 
subleading in l/N^. Both approximations drop only finite terms (in the 1/e expansion) and 
retain the full colour dependence of the IR-divergent pieces. Our explicit setup for the finite 
parts is given in the following. Initially, our virtual contributions are computed in the four- 
dimensional helicity (FDH) scheme 68|, |69|], which is most convenient for the use of helicity 
amplitudes. Only in the end are the virtual amplitudes converted to the more standard 't 
Hooft-Veltman scheme used as input for SHERPA. This conversion is done through 
a simple shift proportional to the born matrix element (see e.g. section six in jl9|). Our 
leading-colour approximation contains the full scheme shift contributions, including terms 
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FIG. 9: A comparison of the full and leading-colour virtual contributions to the pT distribution 
of the fourth jet in VF~'+4-jet production at the 7 TeV LHC. Their difference (full-minus- leading- 
colour) is also shown. We show not a physical cross section, but rather a particular piece of it. The 
subleading-colour contribution is suppressed by a factor of approximately uniformly. The 

vertical lines in the center of each bin, indicate integration errors. 

subleading in l/Nc- In both approaches, the partial amplitudes include the leading terms 
in the formal limit Nc ^ 1 and nj ^ 1 with N^/uf fixed. Here, all interference terms of 



there, in 



the non-zero partial amplitudes are kept. To point out one difference to ref. 
a more conventional approach, only the leading-colour interference terms were kept. 

The comparison of the full-colour and leading-colour virtual contributions to the px dis- 
tribution of the fourth jet in iy~+4-jet production at the LHC is shown in fig. [91 Also 
displayed are the subleading-colour contributions by themselves, labelled as "full-minus- 
leading-colour" . In order to obtain the full parton level differential cross-section one must 
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add the real and born contributions in the usual way. The complete pr distribution, includ- 
ing born, real and leading-colour virtual contributions, has been given already in 

We confirm that the subleading-colour contribution is suppressed uniformly over the 
range from 20 GeV to 200 GeV. The suppression appears consistent with the expected factor 
of l/N^ with Nn = 3. With the leading-colour virtual part accounting for about 20% of the 

nn 

leading-colour total cross section, as already observed in refs. [5|, HlJ , the subleading-colour 
virtual part amounts to less than a 3% correction. Although numerically similar, the present 
leading colour approximation and the one from ref. ll(| are not identical. We have checked 
that our conclusions hold for both definitions, namely that sub-leading colour contributions 
are at the level of a few percent. 

A possible exception to the uniform suppression are zeros of the leading-colour virtual 
cross section. Such zeros are not excluded on general grounds, as the virtual piece alone 
is not physical. A priori, there is no reason to assume that the vanishing of the leading- 
colour contribution forces also the vanishing of the subleading-colour contribution. Close to 
a zero of the leading-colour contribution we naturally expect a relative enhancement of the 
subleading-colour piece. In the range of the distribution presented in fig. [9]we do not observe 
such behaviour and observe a uniform suppression of the subleading-colour contribution. We 
leave a more thorough consideration of this issue, as well as a more detailed phenomenological 
analysis, to future work. 



V. CONCLUSION 



Recent years have seen impressive progress in NLO parton level computations relevant 
for LHC physics. In particular, the unitarity and on-shell approaches have opened new 
possibilities for the computation of one-loop scattering amplitudes. 

An important ingredient in these advances is a clear theoretical understanding of the 
amplitudes, which enables an efficient structuring of the calculation. This progress has made 
possible high multiplicity parton-level predictions at NLO for generic physics processes and 
models. The colour organisation discussed here is one important aspect towards automation 
in the colour-ordered approach. 

We have presented explicit expressions for the colour decomposition of pure-QCD am- 
plitudes with up to seven partons. We also gave their extension to include a leptonically 
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decaying ly-boson. Whereas leading-colour expressions are simple to write down, terms sub- 
leading in 1/Nc are not easily obtained. We have described a general algorithm for obtaining 
such subleading terms, applicable to cases with any number of fermion pairs. 

Although we give explicit expressions, typically these are not unique. There exist relations 
between primitive amplitudes which originate in the anti-symmetry of the fermion-fermion- 
gluon colour-ordered three-point vertex. We expect that these "fermion-flip" relations will 
be important in obtaining all-ra formulae for multi-quark amplitudes, as exist already for 
the pure-gluon 17 1 and two-quark [1^ cases (see also {gOI). 

Contributions to the virtual corrections at subleading order in llN.. have been shown to 

h 

be uniformly small over phase space for W+ 1, 2, 3-jet production in [5|]. However, no general 
theorems are available to assess the size of subleading-colour terms in all regions of phase 
space a priori. We have applied our new results to understand the effect of subleading- 
colour virtual contributions to distributions oi W + 4-jet production. We showed that 
for the fourth-jet pt distribution the effect is small (a few percent) and observe a similar 
behaviour for other distributions. A detailed study of the physics, including subleading- 
colour contributions, to pure jets and vector-boson plus jets could be carried out in the 
future using the results of this article. 

In summary, we have presented a general method towards a colour-ordered approach 
for loop computations with many coloured final states. We expect these results will aid in 
physics studies of relevance to the LHC in the near future. 
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